Is there a Glass Transition in Planar Vortex Systems? 
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The criteria for the existence of a glass transition in a planar vortex array with quenched dis- 
order are studied. Applying a replica Bethe ansatz, we obtain for self-avoiding vortices the exact 
quenched average free energy and effective stiffness which is found to be in excellent agreement with 
recent numerical results for the related random bond dimer model Q. Including a repulsive vortex 
interaction and a finite vortex persistence length £, we find that for £ — > the system is at all 
temperatures in a glassy phase; a glass transition exists only for finite £. Our results indicate that 
planar vortex arrays in superconducting films are glassy at presumably all temperatures. 



The existence of a vortex glass (VG) in type-II super- 
conductors is of crucial importance for their performance. 
Quenched disorder affects both the translational order 
of the Abrikosov vortex lattice and the low current re- 
sponse 0] . The thermodynamic response of a planar VG 
in a mesoscopic superconducting film || has been stud- 
ied in a recent experiment 0], providing 'fingerprints' of 
the sample's disorder M . Glassy phases are expected to 
be generic to a wide class of disorder dominated systems. 
The planar VG is a famous starting point within this class 
since it is among the very few systems where analytical 
progress seems possible. The situation is here analogous 
to that of exactly solvable ID quantum systems which 
have advanced the understanding of strongly correlated 
systems. 

Thermal fluctuations compete against the VG, raising 
the question of the existence of a glass transition at a 
finite temperature T g However ,thc literature on this 
subject is rather contradictory jj], 0, f§. Major support, 
both analytical and numerical, for a finite T g stems from 
the phenomenologically related random field (RF) XY 
model which has a finite T g transition 0, [u^. However, 
considering the entropic vortex interaction, the absence 
of a finite T g transition has been predicted under cer- 
tain assumptions as well The theory underlying the 
vortex system can be applied to a broad class of other 
systems like surface steps jy] , polymers [ 12 and domain 
walls fll3|| . A more fundamental aspect of the planar 
VG is related to its often debated replica description, 
including the concept of replica symmetry breaking, and 
the non-trivial exchange of replica and thermodynamic 
limit. Due to the relation to the recently simulated ran- 
dom bond dimer model 0, the VG provides a unique 
system to test replica theory and could be useful in clar- 
ifying important concepts in the statistical mechanics of 
disordered systems. 

In this paper, we analyze the conditions under which 
a planar vortex lattice undergoes a finite T g glass transi- 
tion. For vanishing vortex persistence length £, we find 
the lattice to be glassy at all temperatures with an an- 
alytic free energy. Analyticity of the free energy is also 
obtained for the vortex-free RFXY model. Our replica 
symmetric Bethe ansatz results are found to be in excel- 



lent agreement with recent simulations Q. For finite £, 
there exists a finite T g which, however, is presumably too 
large to be observable in superconducting films. 

Vortex system. — We consider a lattice of self-avoiding 
directed elastic lines of density p — I/a in 1 + 1 dimen- 
sions. Each line is characterized by its position Xj(z) 
and line tension g. The lines interact via a short-ranged 
repulsive pair potential U(x). Quenched disorder cou- 
ples locally to the lines via the random potential V(r), 
r = (x, z), which is assumed to have short-ranged corre- 
lations V(r)V(r') — A(5(r — r'). The total energy reads 



H = dz 




■ Xj) + V(xj,z) 
(I) 

This model can be considered as the continuum version 
of lines placed on the bonds of a rectangular lattice [|l4j . 
Below, we will also consider lines with a finite persistence 
length £, corresponding to correlated random walks in the 
lattice model. 

Random field XY model. — The line lattice is usually 
treated within a 2D elastic theory for the displacement 
field u(r) so that the line positions are Xj(z) = ja + 
u(ja,z). The effective Hamiltonian is then 



H 



XY 



d 2 r{ C f(d x uf + C f(d z uf + p(r)V(r)} (2) 



with compression modulus en = all" (a), tilt modulus 
C44 = g/a and line density p(r) = ^ - 8{x — Xj(z)). 
Rewriting p(r) as a Fourier series in u(r), the model be- 
comes equivalent to the random field XY model without 
vortices Q. It is well known that, for weak disorder, this 
model shows a phase transition to a glassy phase at a fi- 
nite temperature T = T g . A renormalization group (RG) 
analysis yields the scaling dimension of disorder 



A A = 2 1 - 
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(3) 



giving T g = a JcucaII -n . An important property of the 
XY model is that the elastic moduli remain unrenormal- 
ized due to a statistical symmetry M. Thus, naively 



2 



one would expect for the line lattice a transition at a 
T g determined by the elastic moduli on the length scale 
where Hxy becomes applicable. However, as we will 
show below, there is, in fact, a renormalization of cu in 
the original model of Eq. (Q) which can even destroy the 
transition. 

What are the manifestations of the glass transition in 
the XY model? Anomalous variations in the response to 
an external tilt of the lines have been used to characterize 
the glassy phase [Q . However, signatures of the transition 
in the free energy itself have not been studied so far. 
Using the known mapping of the XY model onto a vector 
Coulomb gas one can derive the RG flow of the free 
energy as shown in for the pure XY model. From this 
analysis we find that the quenched average free energy is 
analytic at T g . Weak disorder makes the contribution 




A 



c?/ 2 c^V(1-t)2(1-2t) 



(4) 



to the free energy density for small r = 1 — T/T g , where 
we have neglected T-independent terms. 

Another signature of the glassy phase are its correla- 
tions §. For T > T g , one finds [((u(r) - u(0)) 2 )] ~ 
log(|r|/a), and in the glassy phase for small r > the 
slightly faster growth [((u(r) — u(0)) 2 )] ~ log 2 (|r|/a). 
Here (. . .) and [. . .] denote thermal and disorder average, 
respectively. The decay of the corresponding density- 
density correlations suggests that the correlation length 
is infinite on both sides of the transition which is consis- 
tent with the predicted analyticity of free energy. Next 
we will focus on the correlations of thermal fluctuations 
Su(r) = u(r) — (u(r)) around the state pinned by disor- 
der. On both sides of the transition one has 



C(r) = [((Su(r) - Su{0)) 2 }} 



lQirK 



log(|r|/a), (5) 



with K = (n/16)T g /T = (a 2 /16)^c llC44 /T. Below, we 
will calculate K directly in terms of the parameters of 
the original model of Eq. (|l|). 

Self-avoiding vortices. — A generic but simpler model 
is obtained if the only vortex interaction of Eq. (|l|) is 
the one representing the non-crossing condition, i.e., a 
hardcore repulsion U{x) — cS(x) with c — * oo. Using 
the replica method, with n replicas, the model can be 
mapped to SU(n) fermions and thus can be solved ex- 
actly by a replica Bethe ansatz (RBA) without replica- 
symmetry breaking []l4|| . To date, only deep inside the 
glassy phase, e.g., for H — > H c i, the quenched average 
free energy density [/] has been studied by this technique 
fl4[ [T^ |. However, at arbitrary H (or density p) the RBA 
yields 



rr2 pQ 



(q)dq , 



(6) 



where /o is the single vortex free energy per unit length; 
its precise form is unimportant in what follows. g(q) is 
determined by 
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Mq-q') 
with length scale 



7T coth (irl d (q - q')) } g{q')dq' 



(7) 



d _ T 3 /gA, and Q is fixed by 
J-q £?(<?) ^5 = P- The integral equation (]?]) can be solved 
perturbatively in Ql d , leading to 



Q(q) = Vi - (q/Q) 2 
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Higher order terms can be easily calculated and will de- 
pend also on q, thus changing the functional form of the 
lowest order result. In the opposite limit, Qld — » oo, 
i.e, for vanishing disorder, the solution is g(q) = l/(27r). 
With the result for g(q) at hand, both pld and the in- 
tegral of Eq. (||) can be obtained as a series in {Qld) 2 ■ 
However, we are actually interested in eliminating Q in 
order to calculate [/] in terms of p and Id- At this stage, 
a particularly interesting property of g(q) is important. 
The exact result for the integral of Eq. (^J) as a function 
of pl d is obtained from the two lowest order terms of Eq. 
(||). Indeed, it can be shown order by order in Ql d that 
the exact result is j^q 2 g{q)dk = ^-p 3 + p 2 ^ 1 . This 
result we could also confirm numerically with high preci- 
sion. Thus, the RBA gives for the quenched average free 
energy density the exact result 



\f\=foP+~ P 3 + ^P 2 - 
6 g 2T 



(9) 



For the interaction (last two) terms it is possible to 
give a simple explanation. The first is just the known 
Pokrovski-Talapov interaction from entropy loss due 
to collisions of thermally fluctuating lines in pure sys- 
tems. The second term can be explained by a similar 
argument which gives the interaction in the dilute limit 
p ^ I . It arises from an increase in energy due to col- 
lisions which prevent the lines from gaining their optimal 
energy from the random potential |l4j . It is rather unex- 
pected that the total effective interaction is just the sum 
of these two contributions. As a consequence of that, the 
free energy is an analytic function. Based on our result 
for the XY model, this does not exclude the possibility of 
a finitc-T glass transition. However, its existence can be 
excludedby computing the scaling dimension Aa, Eq. (|^), 
in terms of the macroscopic compression modulus cJx< It 
is given by the effective interaction, = aU" s (a), with 
Ueg/a corresponding to the two last terms of Eq. (||). 
Using C44 — g/a, and 
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we find Aa > at all finite temperatures, i.e., T g = 00. 
With the exact available, we are able to relate the 
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phenomenological stiffness K of Eq. 
scopic parameters of the model of Eq. 



K = 



16 



agA 

TT 2 T 3 
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|) to the micro- 
|), leading to 

(11) 



For vanishing disorder, A = 0, the previously known 
result K = tt/16 is recovered ]T^ . 

There are a few comments in order to clarify the ap- 
parent discrepancy between the exact RBA results for 
self-avoiding lines and the RG predictions for the random 
field XY model of Eq. (Q). To begin with, on microscopic 
length scales there is no direct link between the two mod- 
els since en = due to self-avoidance as only interaction. 
However, the XY model is generally believed to describe 
the vortex lattice at length scales larger than a where a 
finite fluctuation induced Cn exists || That this point 
of view is incorrect for finding out about the existence of a 
glass transition in the vortex system shows the following 
argument. On large but finite length scales L ^> a 2 g/T 
the effective en for A — > can be obtained from the 
analogy of the self-avoiding vortex lines to ID fermions 
flit . Since a finite vortex length L translates into a finite 
temperature for the fermions, the Sommerfeld expansion 
can be used to get the correction Sen — —ag/3L 2 to the 
result of Eq. (|l(]). Thus, on any sufficiently large length 
scale, disorder seems to be irrelevant since then Aa < 0. 
In contrast, the RBA shows that disorder is always rele- 
vant. The contradicting results are explained by the fact 
that the statistical symmetry in the XY model is broken 
on all length scales by the non-crossing condition. 

Dimer model. - The self-avoiding vortices can be 
simulated directly, i.e., not via the XY model, using a 
mapping to a discrete dimer model with quenched bond 
disorder Q. A recently developed polynomial algo- 
rithm allows to study large systems at finite temperature 
but fixed density p — 1/2 (in units of the dimer lat- 
tice constant) Comparison of our RBA predictions 
with the results of this numerical approach provides a 
critical test of the replica approach, including replica- 
symmetry and proper analytical continuation while tak- 
ing the number of replicas to zero. The dimer model 
is defined as follows: Choose a subset (whose elements 
are called dimers) of the bonds on a square lattice, such 
that every lattice site is touched by exactly one of these 
dimers. The energy of one such dimer covering D is de- 
fined by Hd = J2u j)eD e ij/Td, where the sum is over all 
dimers of D, are random bond energies drawn from 
a Gaussian distribution with zero mean and a variance 
of one, and Td is the dimer temperature. Via the well- 
known mapping to the solid-on-solid model, it is possible 
to assign every dimer covering a height function which 
maps to the line displacement u(r) in the vortex sys- 
tem representation. This representation results from the 
superposition of an arbitrary covering D and a fixed ref- 
erence covering that specifies the direction of the lines, 



cf. Ref. [jD. By measuring the height configurations, Zeng 
et al. could determine the correlation function of Eq. 
and thus the stiffness K as a function of Td H . In order 
to compare their results to our prediction of Eq. ( |Tl] ) , the 
parameters g/T and A/T 2 have to be determined for the 
dimer representation of the vortex system. It is easy to 
see that the ratio of horizontal and vertical line steps on 
the lattice yields g/T — 2. The precise relation between 
A and Td can be obtained by comparing the annealed av- 
erage free energy of the continuum model of Eq. (|lj) and 
the discrete dimer version of the vortex lattice. Using 
known methods for computing the partition function of 
dimer models [[t9) , we get 
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arctan(x) 



(12) 



with £d a regularization length for the disorder correla- 
tor (2(J, and Catalan's constant G = 0.915966. With 
the above relations and a = l/p = 2in Eq. (|ll]) , we ob- 
tain K(Td) as plotted in Fig. |l|, together with the data 
provided by the dimer algorithm of Ref. Q . Obviously, 
there is excellent agreement between our RBA results 
and simulations. 




FIG. 1: Replica Bethe ansatz result for the stiffness K, 
Eq. ([ill), plotted as function of the dimer temperature Td, 
using Eq. ( fl^ ) with £d = 1.15, and the corresponding data 
(circles) from the dimer model algorithm of Ref. . 

Finite interaction and persistence length. - So far, 
we have shown that the model of Eq. (P with a vanish- 
ing single line persistence length £ and U(x) represent- 
ing only a non-crossing condition is always in the glassy 
phase. We now consider the effects of an arbitrary re- 
pulsive interaction U(x) and a finite persistence length 
£. For the vortex system, U{x) diverges for x — ► and 
decays exponentially for large x; £ can be identified with 
the coherence length of the superconductor. 

Let us first discuss a general interaction while keeping 
£ = 0. The relevance of disorder now depends on the 
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interaction induced change in Cn for A -> 0, cf. Eq. (§). 
This change can be obtained by employing the known 
mapping [jl3| of elastic lines to ID spin-less fermions 
which however interact via U(x), denoting now the fi- 
nite distance interaction on top of the self-avoidance. 
The compressibility n of the fermions is related to the 
compression modulus of the vortex system according to 
1/k = a 2 c\\. For weak interaction, k is easily calcu- 
lated, e.g., by bosonization. Using the standard notation 
of Luttinger liquid theory pjj |, one gets K = 1/ttHvn 
with v N = v F + (U(0) + U(2n/a)/2)/hn, where U(q) 
is the Fourier transform of U(x). Thus, any repulsive 
interaction increases the effective compression modulus 



Cxi 



with 



1 



1 - 



-3 U(0) + -U(2n/a) 



(13) 



and for U(x) — > 0, the known Pokrovsky-Talapov result 
cn = c^j T = n 2 T 2 j ga? is recovered. Thus according to 
Eq. (^), the scaling dimension of disorder is increased 
by any finite vortex interaction, and the system remains 
always in the glassy phase. 

Next, we neglect interactions for the time being and 
now consider the thermally fluctuating vortex lines as 
correlated random walks, i.e., we assume a finite persis- 
tence length £. These correlations increase the fluctua- 
tions of a single line. Denoting the so decreased macro- 
scopic line tension again by g 7 entropy reduction due to 
line collisions leads not only to the Pokrovsky-Talapov 
result but also to an additional renormalization of the 
modulus cn = — where, for £ <C a 2 g/T, 
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For larger £, cn increases relative to this first order result. 
Including the effect of disorder on cn, cf. Eq. ([To|), one 
obtains from Eq. (|^) the finite transition temperature 



T g = 0.42 



2 A 



P 3 £ 
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(15) 



Thus, a finite persistence length yields a finite tempera- 
ture glass transition. 

Finally, we study the combined effect of interactions 
and a finite persistence length. If disorder is strong com- 
pared to the interaction, A > T(U(Q) + U(2tt/o,)/2), the 
result of Eq. ( |l5|) remains applicable. In the opposite 
limit, 



Tg is independent of disorder, 



T a = 0.31 



g 2 (U(0) + U(2ir/a)/2) 
P 3 C 



1/3 



(16) 



For vanishing persistence length, T g diverges, and the 
system is glassy at all temperatures. 

Discussion. — The model of Eq. (pi) applies to vor- 
tex systems in films of thickness d <~ a, A with A the 



penetration depth. Then the crucial questions is if the 
glass transition temperature T g can be smaller than the 
critical temperature T c of the superconductor. Since 
the only mechanism leading to a finite T g is a finite 
persistence length, we maximize the effect from in 
Eq. (|l4|) by considering small a <C A. Then Ginzburg- 
Landau theory yields a U (<?) Jp) which is too large for the 
bosonization result of Eq. ( pL3[ ) to be valid. However, a 
lower bound is cj]* > bc\i with 6 of order one. Using 
Eq. fldj), we find T g > 36/10 x a 2 g/£. With the stiffness 
g = <f>o/167r x d/a\ 2 , $ = hc/2e, ||, we obtain in the 
limit of infinitesimal disorder T g > $q/167t 2 x ad/£\ 2 
which becomes, with a, d > £, T g /T c > 1 — 5 • 10~ 4 for 
T c = 100K and with the length £ = lOA, A = 100A at 
T = 0. Thus, T g is extremely close to T c , and presum- 
ably any small amount of quenched disorder will raise T g 
above T c , making it impossible to observe a glass transi- 
tion in planar vortex systems. 
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